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Abstract 

In this work we present a very simple approach to input-output relations in 
optical cavities, limiting ourselves to one- and two-photon states of the field. After 
field quantization, we derive the non-unitary transformation between Inside and 
Outside annihilation and creation operators. Then we express the most general two- 
photon state generated by Inside creation operators, through base states generated 
by Outside creation operators. After renormalization of coefficients of inside two- 
photon state, we calculate the outside photon-number probability distribution in a 
general case. Finally we treat with some detail the single mode and symmetrical 
cavity case. 



1 Introduction 



The problem of interaction between a pair of atoms in free space and in a cavity has 
been subject of several investigations in the past years [|I|, ^, |3|, f|, [7|, §]. Recently 
the spontaneous emission of a pair of two identical atoms or molecules in a planar 
Fabry-Perot microcavity has been subject of theoretical and experimental research 
H [ll]]. This work starts from our attempt to give a simple interpretation to some 
recent experimental results. Consider a process of spontaneous emission of a pair of 
photons by two distinct molecules inside a microcavity. If we measure the number 
of photons emitted outside the cavity, and if there are not dissipative phenomena, 
we will find only three possible results: two photons detected on the right and no 
one on the left, two photons detected on the left and no one on the right, one photon 
detected on the left and one on the right. In this work we look for a solution to this 
question: if we know the distribution of the number of photons outside the cavity, 
how can we obtain information about the process of emission which generated the 
photons without completely solving the problem? Our basic idea is to describe the 
electromagnetic field inside the cavity by means of a two-photon state as general as 
possible; the coefficients of the expansion of this state in a proper basis set depend 
on the process which generated the state itself. If we project this state on the basis 
number-states defined outside the cavity, we obtain directly the probability distri- 
bution we look for. Therefore if we change the above mentioned coefficients, we can 
study how the probability distribution changes, and by comparison with the mea- 
sured distribution we can obtain the values of these coefficients getting information 
about the process active inside the cavity. This procedure is not orthodox or fully 
justified, however the obtained results are in qualitative agreement with experimen- 
tal results. 

The relations of Input-Output and their connection with the traditional stochastic 
methods based on Langevin equations have been studied by Knoll et al. fl22] . The 
aim of the present paper is to generalize the methods used to describe an optical ele- 
ment with two inputs and two outputs, as a Beam-Splitter, to derive simple relations 
between fields inside and outside a planar Fabry-Perot cavity. In Sec. II we quantize 
the electromagnetic field generalizing to three dimensions an approach presented by 



Barnett et al. 11 for one-dimensional fields. In Sec. Ill we restrict our attention 



to one-dimensional fields and derive the relations between operators defined in the 
space Inside and Outside the cavity. Then, in Sec. IV we build and study the 
states of the field generated by linear and bilinear forms of creation operators both 
inside and outside the cavity and, in Sec. V, we calculate the outside probability dis- 
tributions for these two-photon states. Finally, we summarize our results in Sec. VI. 
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2 Spatial modes and Field Quantization 



Now we calculate the appropriate normal modes for quantization of the electro- 
magnetic field. In a traditional approach, one first determines the modes of the 
classical boundary value problem, then one quantizes the field in terms of these 
modes [O, fTl. An alternative approach has been presented by Barnett et al. [TT 



We have generalized this work, restricted to an one-dimensional fields, to three- 
dimensional fields. 

In free space, having fixed the Coulomb gauge, the spatial variations of the electro- 
magnetic field can be described by the solutions Ufc(r) of the Helmholtz equations 



W(r) + fc 2 u fe (r) = 0. (1) 

For each value of k, there are two solutions that we indicate with iiR A (k, r) and 
UL\(k., r). We choose them as plane waves that are incident respectively from left 
and right on the plane of equation z = const., as shown in Fig. 1. 
Here A = 1, 2 is the polarization index and k is the wave- vector such as |k| 2 = k 2 = 
u 2 /c 2 . The explicit form of UR A (k, r) and UL A (k, r) is therefore 

URx{k, r) = e A (k + ) exp(ik + • r), 

(2) 

u LA (k, r) = e A (k_) exp(ik_ • r), 

where we have defined 

k± = fc(sin# cos0, sin 8 sin0, icos^), (3) 



for (0 < 9 < tt/2), and 

6i(k±) = (sin0, — cos0, 0), 



(4) 

62 (k±) = (cos^ cos0, cos^ sin <fi, =)= sin 6). 
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FIG.l. Scheme of the plane waves u A (k, r) and u' A (k, r) incident respectively 
from left and from right on the plane z — 0. 

Now, following a standard procedure |15|, [Uj, we introduce the mode creation and 
destruction operators a# A (k) and OL A (k) as 



d^, A (k) and d_R A (k) for the mode UR A (k, r), 
d^ A (k) and flLA(k) for the mode ii£ A (k, r). 



(5) 



Being k a three-dimensional continuous variable, these operators satisfy the follow- 
ing commutation relations: 



(6) 



d RA (k), aJ ZA ,(k') = d LA (k), a\ x ,Qs!) = <5 AA /<5(k - k'), 
a LA (k), dJj A ,(k') = a/j A (k), a[ x ,(k') = 0. 
The vector potential operator is written, in Heisenberg representation, as 

A(r,t) = A + (r,t) + A-(r,t), (7) 



where 
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X 



^2 \ u rx(K r)a RX (k) + u LA (k, r)a LA (k) exp(-iwt) 



A=1.2 



and A"(r,t) = [A+(r,*)]t. 

Now we consider an infinitesimally thin non-absorbing dielectric slab, placed in the 
plane z = zq and represented by a non-homogeneous dielectric constant of equation 



e(z) = e [t]5(z - z ) + 1] 



(9) 



where r\ is a constant depending upon the optical properties of dielectric media (17 . 



a L \(k) 




b R x(k) 




FIG. 2. Schematic representation of the dielectric slab that transform Input 
operators a R \(k) and l>L\(k) in Output operators a L x(k) and b R \(k). 



The complex reflection and transmission coefficients r\(k) and t\(k) of the slab 
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depend upon the incident radiation wave-vector and polarization and they have the 
following unitary lossless properties, valid for all values of k, A: 



|r A (k)| 2 +|t A (k)| 2 = l, 

(10) 

t A (k)r A (k)+r A (k)t A (k)=0. 

Consider the arrangement in Fig. 2, which shows the interaction of two fields of 
wavevectors k + and k_ respectively, on a dielectric slab placed at z — 0. It is well 



known [L9[] that the annihilation operators a_L A (k) and CR A (k) associated to Out- 
put modes are related to the annihilation operators &R A (k) and c^ A (k) associated to 
Input modes, by the 2x2 matrix 



'a LA (k)\ /r A (k) t A (k)\ /a m (k) N 
v Siix(k)J U(k) r A (k) / Wk), 



whose unitarity is assured by Eqs. fllP)) . Now imagine to put a pair of infinitesimally 
thin dielectric slabs at z = ±1/2, to form a planar Fabry-Perot cavity, as shown in 
Fig. 3. 

We impose boundary conditions of the form ( |iT|) both on slab 1 and on slab 2 ob- 
taining respectively, 



Mk)e^H/ 2 ) = t 1A (k)a /?A (k)e^(- / / 2 ) +r 1A (k)6 LA (k)e^(- i / 2 ), 
a LX (k)e^- 1 ^ = t 1A (k)6 LA (k)e^(-^ + r 1A (k)a RA (k)e^/ 2 \ 



(12) 



and 



c RX (k)e^/» = t 2X (k)b RX (k)e^/V + r 2A (k)c iA (k)e^/ 2 \ 
b LX (k)e^/V = t 2A (k)c iA (k)e^ 2 ) + r 2A (k)6 RA (k)e^/ 2 \ 



(13) 



where ^Pr{z) and ^l{z) are the phases generated by field propagation inside the 
cavity. 



Mirror 1 Mirror 2 




z = -l/2 z = l/2 



FIG. 3. Schematic representation of a planar Fabry-Perot cavity with notation 
for Input operators a RX (k), C£>(k), Inside operators bn\(k), b^k) and Outside op- 
erators a L \(k), cja(k). 

We freely choose the zero phases in the middle of the cavity, so that tpp(-z) = 
—<f F (z), where F = R, L and ^Pr{z) = —<pi,(z). Then we can put 

VL (±l/2) = T 6/2, p R (±l/2) = ±5/2, (14) 

where 

5 = -/cosfl, (15) 

c 

is half of the phase gained in a double traversal of the cavity |[20|| . From the first of 
Eqs. (|12|) and the second of (|T3"D , we can express the intracavity operators bR\(k) 
and bi\(k) as 



6 



b RX (k) = a RX {k)I x {k) + c LX {k)J' x {k), 



hx(V = c LX (k) J A (k) + a RX (k)I' x (k). 



(16) 
(17) 



where we have defined, following Ref. [13 



D x {k) = 1 -ri A (k)r 2 A(k)exp(2z<5), 

I x (k)=t lx {k)/D x (k), 
J A (k) = t 2A (k)/ J D A (k), 

J A (k)=r 2A (k)exp(^)/ A (k), 
J A (k) = r 1A (k)exp(z5)/ A (k). 



(19) 



(20) 



Substituting Eqs. 
obtain 



BP- ( PI) in the second of the (|12"D and in the first of (|13|), we 



c RA (k) = a m (k)T A (k) + a LA (k)i? A (k), 



a LA (k) = c LA (k)i? A (k) + a^k^k), 



(21) 
(22) 



where 



T A (k) = T A (k) 



t 1A (k)t 2A (k) 
Dx(k) ■ 



Rx(k) 



r 1A (k) exp(-i5) + r 2A (k) [t 2 ix (k) - rj x (k)] expjiS) 

Dx(k) 

r 2A (k) exp(-^) + r 1A (k) [t 2 2X (k) - r| A (k)] exp(^) 



(23) 



(24) 



(25) 



The coefficients R\(k) and T A (k) and the corresponding primed represent the reflec- 
tion and transmission coefficients of a Fabry-Perot cavity as a whole. It can readily 
be shown that they satisfy the conditions |13| 



\R x (k)\ = \R' x (k)\ , 



|i? A (k)| 2 + |T A (k)| 2 = K(k)| 2 + |r{(k)| 2 = l, 



i? A (k)7l(k) + # A (k)T A *(k) = 0. 



(26) 
(27) 
(28) 



Using Eqs. (H) and (^)-(^) it is not difficult to show that the Output operators 
CR A (k) and a^ A (k) satisfy canonical commutation rules 



a LA (k), 4 A ,(k') = ^(k), c^(k') = 5 xx ,5(k-k') 



t nj\ 



(29) 



a LA (k), 4v( k ') = c RX (k), a\ w (k 



0. 



while the intracavity operators &R A (k) and &L A (k) satisfy anomalous commutation 
rules 0, [TT 



b LX (k), b j RX ,(k') =S X yS(k-k') 



S X y5(k 


-k') 


1 - 


-|r 1A (k)r 2A (k)| 2 


|1 - 


r 1A (k)r 2A (k)e M | 2 

(30) 


~hx(k) 


&L 


(k')" 




r 2X (k)e u 


[1- 


|r 1A (k)| 2 ]+r* A (k)e- 5 [l-|r 2A (k)| 2 ] 






|1 - 


r 1A (k)r 2A (k)e 2 ^| 2 



h RX (k), b[ y (k>) 



(31) 



These equations are the three-dimensional generalization of Eqs. (9) given in Ref. 
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11] for one- dimensional fields. 

Because of the presence of the cavity, the vector potential is now written as 



A+(r,t) = dk 



h 



167T 3 eo^ 



1/2 



^2 F A (k,r)exp(-z^t), 



(32) 



A=l,2 



where 



ik_ r 



a iJA (k)e A (k + )e jk +' r + a LA (k)e A (k_)e 



F A (k,r) = ^ 6 i?A (k)e A (k + )e jk +- r + 6 iA (k)e A (k_)e jk - r 



[c i?A (k)e A (k + )e ik +' r + c LA (k)e A (k_)e 



ik— r 



-oo < z < 1/2 



-l/2<z<l/2, (33) 



1/2 < z < oo . 



We note that using different pairs of annihilation operators for each region of space 
delimited by the cavity, as in Eqs. ( |3~2"l - j33| ) , we obtain a free-field like representation 
for all space inside and outside the cavity, but the canonical commutation rules are 
lost for intracavity operators, as shown by Eqs. (|30| - |3T]) . Conversely if we choose the 
mode functions for example as in Ref. |14| |, we lost free-field like representation, but 
we obtain canonical commutation rules for annihilation and creation operators in 
whole space. It is easy to show that our result agrees with that of Ref. ]T3| , indeed 
defining 



ora(1<0 = a k A, 
c LA (k) = a' kA , 



(34) 



and substituting Eqs. (TT6I - ITTD and ( f2T[ - |22|) into Eqs. (|33D , by a straightforward cal- 
culation, we obtain 



F A (k, r) = e(k, A) (U kX a kX + U kX a kX ) , -oo < z < +oo, 



(35) 



where now the mode function U k \ are defined differently on the three regions of the 
space, as shown in the tables (2.6) and (2.7) of Ref. O . 
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3 One-dimensional formulation 



In this work we are interested to derive Input-Output relations for a single transverse 
mode of the cavity, having finite cross-section area A orthogonal to the z axis which, 
in fact, depend upon the geometrical and transmitting properties of the cavity itself 
|5|, E|. Then, following Ref. PT[, we impose periodic boundary conditions on both 
directions x and y, so that the corresponding components of the wave-vector are 
restricted to the discrete values 



2vr 2tt 

~JpfL X , ky = 



n 



where n x ,n y = 0, ±1, ±2, 



(36) 



while k z persists to be continuous positive variable. Then the following conversions 
are required: 



dk 



A 



dk z 



(37) 



5 3 (k-k') 



A 



(27T) 



;S nn >5(k z - k' z ), 



(38) 



a RX (k) ^ (2n/A 1/2 )- 1 d RX (2n/A 1/2 n,k z 
c LX (k) ^ (2n/A 1/2 )- 1 c LX (2n/A 1/2 n,k z 



(39) 



where n = (n x , n y ). Now we fix a linear polarization parallel to the x axis, consider 
only field excitations with n x = n y = and define 



a RX (27c/A 1/2 n, k z ) 
c lx (2ti/A 1/2 n, k z ) 



n=0 



n=0 



c 1/2 d R (uj), 



c 1/2 c L H 



(40) 



where we have defined k z = k = u/c. These operators have the commutators 



a R (uj), a R {u') = c l (uj), c\(J) =8(u-J), 



(41) 



c l (uj), d R (u 



0. 
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Until now, we have limited our attention only to incident operators, however re- 
peating the same procedure for operators 6^>(k), b R \(k) and a£,\(k), CR\(k), we 
define straightforwardly the corresponding one-dimensional operators b R (ui), b^cu) 
and a L (u), c R (u). 

Now we assemble the six operators a R {ui), cl{uj), b R (u>), and ax, (a;), c R (cu) in 

three doublets defining 





(42) 



b(u) 





(43) 





l-l 


f c R {u) \ 









(44) 



Their geometric meaning is illustrated in Fig. 4. 



Mirror 1 



Mirror 2 



Si(u;) 



C2M 



&1M 



62 (w) 



z = -Z/2 



z = Z/2 



Ci(cj) 



a 2 (u;) 



FIG. 4. Schematic representation of a one-dimensional Fabry-Perot cavity with 
notation for Input operators di(u), % = 1,2, Inside operators bi(u>), i — 1,2 and 
Outside operators Cj(u;), i — 1,2. 
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Using Eqs. (|I6|-17) in one- dimensional form, it is not difficult to show with a 
straightforward calculation, that the Inside operators bi(u), % = 1,2 are related 
to Incident operators di(uj), i — 1, 2 by the relations 



(45) 



where the matrix elements Bij{uj) of B(u) are given by: 



B«(u})=ti(u)/D{w), 



2 = 1,2, j = 1,2, (46) 



r,(o;) exp(iujl/c)Bjj(uj), 



where 



D(uj) = 1 — T\(u)r2(u) exp(2iul / c) . 



(47) 



Similarly, the Outside operators Cj(cj), i — 1, 2 are related to the Incident operators 
by 



c(u;) = C(a>)a(a;), 



(48) 



where the matrix elements Cij(u) (i,j = 1, 2) of C(u) are given by: 



f C«(u;) = t x (u)t2(u)/D{u)), 



C i:j (u) 



rj{uj) exp(—iul/c) + Tiiyj) exp[iul/c + 2i axgtj(u)] 



D(u) 



(49) 



Using Eqs. (|10|)-([49|) it is easy to show that the matrix C(u) is unitary: 



C(uj)-C\uj) = C\uj)-C(uj) = 1. 



(50) 



12 



This leads, with Eqs. (f41~l)-(f48|), to the result 



(u;),a}(u/) = c^iv), c}(a/) = 5ijS(uj - J), i,j = 1,2. (51) 



For Inside operators the commutation rules are found to be 



[B{uj)B^J)]..5{uj-uj>) 



Gij(u)5(uj - u'), 



where, for construction, Q{oS) = G^{uj), being Gij(u) = [G and 



(52) 



r2(to>) exp(«u;Z/c)[l — |ri(u;)| 2 ] + r\{uj) exp(— iul/c)[l — \r 2 (Lj)\ 



(53) 



Noting that Det[$(u;)] 7^ for t\(ui) 7^ and £2^) 7^ 0, we can invert it to express 
the relation between Inside and Outside operators as 



where 



c(w) = C{uj)B- 1 {uj)h{uj) 



= M(u)b{Lu), 



(54) 



/ 



t*(u;) 



— p-r exp (-lul/c) 



exp (—iujl/c) 



(55) 



It should be noted that the generic 2x2 matrix that represents a Beam-Splitter 
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must be unitary to preserve the canonical bosonic relations of commutation for 
both Input and Output operators; in our case M.(uj) is not unitary at all. In fact, 



while Output operators Ci(u) satisfy the canonical relations (|51|), Input operators 



bi(uj) satisfy the relations fl52|), that is the so-called anomalous relations of commu- 



tation. Using Eqs. ( |54"| - |55| ) we obtain 

M\uo)M{uo) =g~ 1 (u;), (56) 

or, equivalently, 

M(uj)G{u)M\<jj) = I. (57) 



Because of the non-unitarity of A4(u), the photon-number operator is not con- 
served on a single mode. In fact from Eq. (|5"4]-|5"6|) we obtain 



c\u) ■ c(w) = ■ g'\uj) ■ b{u) + ■ b(w). (5* 



Anyway, because we are working with linear transformations, the most general bi- 
linear form in Inside creation operators b\{u) is still bilinear in Outside creation 
operators c\(u): 



7n44 + 7i2cI4 + 72x44 + 72244 

(59) 

=Pnb{b\ + (3 l2 b[b\ + (3 21 b\b[ + faUfrl 

Because the first of the two forms applied to free space generates a two-photon 
state of the electromagnetic field, the equality ensures the same for the second 
one. Therefore, if it is possible to associate to the most general two-photon state 
generated by Inside operators a state with two photons physically generated inside 
the cavity, then, from the relations between Input and Output operators, we can 
obtain information about the field outside the cavity, that is the actual object of 
measurement. In the following section we study how we can do this. 
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4 States of the Field 



We define the states generated by the linear and bilinear forms of Inside operators as: 

6j(w)|0) = iF^y, In), 

(60) 

St(w)6t(u/)|0) ee ^(00), F^oo')- In), 

where is a label which depends on continuous variable 00 and on discrete vari- 

able i — 1,2. Because of (0) these states are not orthogonal: 

(Fi(u); In^u;'); In) = Gij(oo)5(oo - <J), (61) 

and 

(FiM, Fj(oj 2 ); ln\F k (u 3 ), F,(w 4 ); In) = 

Gik(wi)G J 7(w 2 )£( a; i - w 3 )5(u;2 - w 4 ) (62) 
+ G^O^ikM^O^ - ^3)^1 - w 4 ). 



From Eqs. ( |6l| - |62]) we note that anomalous commutation rules, represented by the 
2x2 hermitian matrix Q(oo), form a metric in the two-dimensional Hilbert space 
generated by Inside operators bi(u). For example if we write the most general one- 
photon state created by Inside operators as 



2 r 

/ duKi(u)\Fi(u);1n), 
i=i J 



(63) 



where KAuj) G C, then its norm is 



1,2 . 

(010) = Yl / dooK* (00)0^(00)^(00) 



(64) 



dwK^w) • Q(oo) ■ K(oo), 



15 



where it is clear the metric-like role of matrix Q(uS). 

As before, we can define the states generated by the linear and bilinear forms of 
Outside operators as 



4(w)|0) = |F,H;Out), 
c\{u)c){u>)\U) = \Fi(u), F^u'); Out), 



(65) 



which, using Eq. (fjj), can be written in terms of Inside operators: 



1,2 



|F 4 («);Out) = ^M^HlF.HsIn), 



fc 

1.2 



^(uO.F^o/JjOut) = ^M^M^(^)I^M,^(^);In), 



(66) 



where My(u;) = [.M]^. Of course these state are orthonormal. 

It is possible to make number-states also for a continuous distribution of modes, 



following Blowet a/.'s method pl| . On this purpose we define two operators Ci{rj) 
as 



Ci{rj) = I dujr)*(uj)ci(uj), i = l,2. 



(67) 



where r]i{uj) are two arbitrary complex functions which satisfy the normalization 
condition. 



duj\rji{uj)\ 2 = 1, i = 1,2. 



(6£ 



It is easy to verify that these new operators obey the following commutation rela- 
tions: 



Oij, 



(69) 
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and therefore they can be used to make the number-states by the usual method, 



\Fnv);Out) = (n\)- 1 / 2 [Cj( V )T\0). (70) 
Using flBgf) you can verify that the states \F™(rj) \ Out) are correctly normalized: 

(^(77); OutlF" 1 ^); Out) = <%<W (71) 



The construction of number-states is more problematic for Inside operators. As 
before we define 



A(0= / du${u)bi(u), £ = 1,2 , (72) 



where £i(ui) are two complex arbitrary functions which can be chosen to satisfy the 
four conditions 



= J dwQ^G^u)^) = r y (0, (73) 

where ry(£) is a given matrix. Suppose to set T^-^) = Sy. Because r(£) is hermi- 
tian by construction, Eq ( |T3"D corresponds to 2©2 conditions, two real and a complex 
one , which the two complex arbitrary functions £i(u)) (£ = 1,2) must satisfy. The 
modules can be determined imposing Tjj(£) = 1. In fact, given an arbitrary function 
£(u;) such as 

J ^|e»| 2 = l, (74) 
if Ga(u) 7^ 0, that is if \ri{uj)r 2 {uj)\ 2 ^ 1, we can write as 



= l [ > * = M, (75) 
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where <pi{uo) is an arbitrary phase, and we can obtain Ta{^) = 1. Phases are still 
arbitrary, but in off-diagonal elements of r(£) there is only the phase difference be- 
tween £i(u;) and ^(^) which is not sufficient, by itself, to satisfy the two requested 
conditions. In fact, you can see that only proper linear combinations of Inside op- 
erators can generate canonical commutation relations. Consider the unitary matrix 
U(oo) that makes Q{oS) diagonal 

U\u) -Q{uj)- U{uj) = V(u), (76) 

where the diagonal matrix T>{uS) has elements Dij{uj) = Xi(cu)8ij, being Xi(u), i = 
1, 2 the two G(u>ys eigenvalues, 

A i H=GnH-(-l) i |Gi 2 (w)|, 2 = 1,2. (77) 

Then if we define the operators di{u) as 

d(w) = U ] {uj)h{uj), (78) 

you can see that they satisfy the following "quasi-canonical" relations: 

[di(u), d](u')] = X^SijSiu - u'). (79) 

If we want to obtain fully canonical relations it is necessary to break the unitariety 
of the relation between operators di(u>) and bi(u>) introducing the matrix £{u) with 
elements = (Aj) -1 / 2 ^- and to define the operators b'^oo) as 

b'(u) =£(u)d(u). (80) 

In fact an easy calculation shows that 
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where, using Eq. (|80D, 



Vfr) = ■=== \e^/%{u) - (-lye-^kiu;)} j = 1,2 , (82) 

and 4>{uj) = arg[Gi2(cu)] is the relative phase of the two components of Q{oS) eigen- 
vectors. 

Of course by means of these operators we could make orthonormal Input number- 
state but they will result to be associated to functions sm(ujz/c) and cos(u>z/c) and, 
as a consequence, the free field-like representation will be lost. 

However it is still possible to write the most general two-photon state generated 
by Inside operators bi(u) as, 

\^ = Ylj du j ^%>V)I^M,^V);in>, (83) 

hi 

where, by construction, the matrix K.(u,u') of elements Kij(u,u'), satisfy 

JC{u,u') = JC T {u',u), (84) 

where "T" indicate transposition. In fact the matrix JC(u,u') is fixed by the emis- 
sion process inside the cavity. 

In a simpler way, a two-photon state generated by Outside operators fl67D , can be 
written as 

|F a fo),F t fa);Oiit) = (2- 1 / 2 )^Cl(r ] )Cl(r ] )\0) 

(85) 

= (2-^)^ du dw / % (a;) % (a; , )|F (a;),F 6 (a;');Out), 



where (a, b = 1,2). In the next section we will show how to express this state by 
Inside states. 
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5 Two-photon states probability distributions 



It is well know that the inverse of photon mean flight time in a planar Fabry-Perot 
cavity is given by p3] 



cl - 


n 




r 2 


»| 


" I 2|n 




?~2 


(") 


|l/2 



7cav = 7 o , : n ' (86) 



Now we consider the spontaneous emission of a pair of two identical atoms or 
molecules within the microcavity [fLOR . Let 7 a tom the single atomic decay rate. In the 
atom-dominate decay regime (that is when 7 cav <C 7 at0 m HU ) , for l/7 a tom "C i <C 
l/7cav the electromagnetic field can be found in a state like to If the matrix 
M. (u) should be unitary, we should calculate easily, as in the quantum theory of a 



lossless beam-splitter |24|, the probability distribution of photon number-states out- 
side the cavity. This is not our case, however we will see that after renormalization 
of the state \ip) coefficients Kij(uj,uj'), it is possible to obtain significant results. For 
this, we calculate the probability amplitude to find the electromagnetic field, rep- 
resented by the state \ip) within the cavity, in the outside state \F a (r]), Fb(r]); Out), 
(a, 6 = 1,2). It is simple to show with the use of Eqs. ( |66] ) and (p2|)-(p5|), that result 
is 



J7) 



(F a (r]),F b (r)y,Out\iP) = 2(2~ 1 / 2 ) 5 * i ' J dw J du'r) a (u)r) b (u;') 

x [M(u) ■ G(u) ■ K{cu, u') ■ Q t (uj') ■ M t (uj')] ab . 
By a lengthy but straightforward calculation, it is simple to show that 

[M(u) ■ G{u) ■ K(u,u/) ■ £V) ■ M T {to')\ ab = Pab(u,u/), 

where we have defined the 2x2 matrix elements P a b(uj,uj') as 
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P u (u,u') 
P 22 (uj,lu') 



£ 2 (u;)£i(u/) [K n + K 2 2a 1 (to)ai(uj') + K 12 a 1 (to') + K 21 o>i(uj)} 
C 2 (uj)Ci(uj') [K n a 2 (uj') + K 22 ai(uj) + K 12 + K 21 ai(uj)a 2 (uj')} 
Ci(uj)C 2 (uj') [K 22 ai(uo') + K n a 2 (uj) + K 21 + Ki 2 a 2 (uo)ai(uj')} 
£i(u;)£ 2 (u/) [K 22 + Kiia 2 (uj)a 2 (uj') + K 2 ia 2 (u') + K 12 a 2 (u)} 



being 



£, 



1,2. 



(90) 



Now we evaluate the ratio TZo ut {R, L\R, R) between the probability Po nt (R,L) of 
observing one photon behind mirror 2 and one photon behind mirror 1 (coincidence), 
and the probability Po u t{R, R) of observing two photon behind mirror 2. 
Using Eqs. ( |57D we obtain 



T^OntyR, L\R,R) 



PputjR, L) 
Pout(R, R) 



(F 1 (r l ), F 2 (r,); Ont\ 



(F?(v);Out\ 



J 


f du>i j 


f du 2 r]l(uJ 1 )r] 2 (uj 2 )P 12 (uj 1 ,uj 2 ) 


2 -v> 


f du. 


'J 


f du 2 rf x {uJijil MPu (<*>i, w 2 ) 



(91) 



We now illustrate the meaning of this formula. We start writing explicitally the 
value of the ratio Pi 2 /Pu\ 



Pi2(vi,U2) = K u r 2 {u; 2 )e iu;2l/C + i^n^Qe^ + K l2 + K 21 r 1 (a; 1 )e^ 1 '/ C r 2 (a; 2 )e^'/ C 

P u (wi,w 2 ) ~ AT n + K 22 r 1 (a; 1 )e^ li / c r 1 (a; 2 )e iW2 '/ c + f^n^e^ + X 21 r 1 (a; 1 )e ia;i '/ C 

(92) 

This expression is only apparently complicated, but each term at numerator and 
denominator is susceptible to a clear physical interpretation. We assume that 
Kij (u>i, uj 2 ) is proportional to the probability amplitude that a pair of excited molecules 
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within the cavity, emit spontaneously one photon with angular frequency lo\ on mode 
i, and one photon with angular frequency 0J2 on mode j, being the proportional- 
ity factor the same for all coefficients K^{oj\, to?). More precisely we assume that 
\Kij(uJi, u>2)\ 2 dujiduj2 is proportional to the emission probability of one photon on 
mode i with angular frequency between u)\ and iO\ + du)\, and one photon on mode 
j with angular frequency between u 2 and u 2 + du 2 - At this point it is easy to see 
how each term which appears in the ratio (^) admits a clear physical interpretation. 



K12 + K2ir 1 (LU 1 )e^ ll/c r 2 (iU2)e iuJ2l/c 



+ 



K u + K 2 2r 1 (iUi)e iuJll / c r 1 (uj2)e 



D iu>2l/c 



+ 



Kun^y^ + K 21 r 1 (uj 1 )e^ l / C 



+ 



FIG. 5. Diagrams illustrating the probability amplitudes (reported in the left 
column), relative to Eq. (^) Here ri(u) [r 2 (u))] is the reflection coefficient of mirror 
1 (at the left) [2 (at the right)]. The photon of angular frequency LO\ is always plot 
higher then photon of angular frequency uo 2 . 
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With the help of Fig. 5 we can see, e.g., that the first term at numerator, corre- 
sponding to the first diagrams of the second row, give us the probability amplitude 
of simultaneous emission of a pair of photons toward right, and that the photon of 
angular frequency uj\ is detected behind the mirror 2, while the photon of angular 
frequency u 2 is detected behind the mirror 1 after reflection on the mirror 2. The 
transmission coefficients and all contributes generated by multiple reflections on the 
cavity mirrors, are computed into terms £2(^1)^1(^2), which we have simplified into 
Eq. (P2]). All the other terms in Eq. (^), admit analogue interpretation shown by 
remaining diagrams in Fig. 5. 

Of course, for reasons of internal consistency of the theory, it need renormalize the 
coefficients K^j (0^1 , uj 2 ) imposing 

1 2 

Y^jduxj duzlK^,^)] 2 = 1. (93) 

At this point it is easy to obtain the correct probability distributions. If we define 

R-Out(R, L\R,R) =TZ\, 

(94) 

T^Out(-R ; L\L,L) = TZ 2 , 
and we impose the normalization condition 

Po nt {R, R) + Pout(R, L) + P 0nt {L, L) = 1, (95) 
the desired distributions are then obtained after some algebra in the form 

Pout(R, R) 



n 2 


n 1 - 






K X K 2 




VR 2 + TL{R 2 







Pont(R, L) = - : ' ' „ , (96) 



Pout(L, L) 



n 1 + n 2 + H{R 2 



where, for example, Po nt {R,R) is the normalized probability to find two photon 
outside the cavity behind the mirror 2 of the cavity. 
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5.1 Single mode 



Now we suppose that the field mode spectrum is discretized by an appropriate pro- 
cedure |?T|, furthermore we fix the attention on a single mode of assigned angular 



frequency to. The commutation relations for creation and annihilation operators 
defined on this discrete set of modes, are written as 



[cii,a]] = [ci,cj] = Sij, (97) 



[kM} = G ir (9* 



Now we can have two photon on a single discrete mode, therefore we define the 
states generated by bilinear and quadratic forms of Inside and Outside operators, as 



and 



\Uf 3 M^(^ 1/2 ) 5lJ m]\0), (99) 



(100) 



|/ i ,/ i ;Out)^(2-^)^clc» 

= £ (2-V»)-'« M^lh, In), 

fe,/ 

respectively. It is easy to see that 

(/«, /,-; In|/ fc , /,; In) = ^-V*) 4 ^ [G ik Gji + G«G iJfc ) . (101) 

Exactly as before, we define the most general two-photon state created by Inside 
operators as 



1,2 

Y^Kylhfoln), (102) 
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where fC = /C T , and we calculate 



(f a , f b - Out|V>> = 2 (2-^) 5ab (M-g-JC-g^M?) 



ab ' 



(103) 



where we have defined the matrix /C as 



(104) 



Therefore Eq. is formally still valid if we make the substitution — > a/2. 
From this point we consider the case of a symmetrical cavity, that is we assume 
ri(u) = r 2 {oj) = r(u) and t\(uj) = t 2 (u) = t(cj). For simplicity we choose the phase 
of transmission and reflection coefficients, as 



t(uj) = iy/l — R and r(uj) = — \/R, 



(105) 



and redefine, into a more expressive form, 



Ku = C RR , 



K 22 = Cll, 



K 12 + K 21 = 2K 12 = C RL . 



(106) 



Then in discrete mode representation, the normalization conditions, can be written 

as 



Crr\ + \Crl\ + \C, 



'LL 



1. 



(107) 



Finally we can write: 



7?Out(-R, L\R,R) 



V2(C RR + C RR )re^ c + C RL (1 + r 2 e 



Crr + C LL r 2 e^ + V2C RL re^ 



:io8) 



We note the presence of factors \/2, which we have introduced because of different 
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normalization required by discrete mode spectrum. Because of them, there is no 
exact correspondence between diagram of Fig. 5 and terms into Eq. ( |108| ). This 
factors arise since we are not using an orthonormal base and therefore a mixing 
between normalization of states with a photon on mode and two photons on mode, 
is produced. Indeed we will see, in the next section, that working with single photon 
states which admit only a single normalization factor, it is possible to obtain a direct 
association between diagrams and formulas. 

Now we calculate Eq. ( [L0S| ) for some particular value of the state \ip). Let co- 
inciding with each of the three states of the orthonormal base defined in Appendix. 
It can readily be shown that 

a) |^) = | n± ) 

C RL = ±^=, C RR = C LL = \ n Qnt (R,L\R,R) = 2. (109) 



b) |V> = |n ) 

C RL = 0, C RR = -C LL = -^= K Ont {R,L\R,R) = 0. (110) 

It is remarkable that for these states having high symmetry, TZont does not depend 
either on mirrors reflectivity, nor on the phase ul/c. Following our interpretative 
scheme, Eq. (|109|) shows that when the emission probability of the pair of photons on 
the same way or on the opposite way, is the same, that is IC^ 2 = \C RR \ 2 + \Cll\ 2 = 
1/2, the probability to observe a coincidence is twice with respect to the probability 
of not observing. Instead when the two photon are emitted along a common way, 
but with pair emission probability amplitude toward right and left which differs for 
a sign, Eq. ( |110| ), the probability of observing a coincidence is zero. But within 



our interpretative scheme, which require none distinction between left and right for 
emission of the pair of the photon in a symmetrical cavity, it is hard to think that 
this state really exist. Therefore a state |no) having C RR ^ Cll is difficult to accept. 
From this point we consider only the case C RR = Cll and we define 



^ = C = IC|e JargC . (HI) 
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Then using Eq. ( |111|) , Eq. (|108|) can be written as 

"^•Out(-R) L\R,R) = 

,— (112) 
8R 2 z 2 - AV2R [cos(x + y) + R cos(x - y)) z + 1 + 2R cos 2x + R 2 



(1 + 2.RCOS2X + R 2 ) z 2 - 2y/2R [costs - y) + i?cos(:r + y)] 2 + 2R' 
where we have used the following notation 

x = ul/c, y = arg(, z=\C\. (113) 



This expression seems still rather complicated, but we can learn something from it 
considering the limit cases Crr = -v4> z = and Crl = 2; = 00: 



, ^ ,~ rln „ N 1 + 2i?cos2a; + J R 2 _ 
limftout = ^5 = ft , 



(114) 



8i? 2 

lim n 0nt (R, L\R,R) = = K c 

z-*oo ° utv ' 1 ' ; 1 + 2.RCOS2X + .R 2 



The first of Eqs. (|114 ) is shown in Fig. 6. From Fig. 6b we can see that when 



R — > 0, TZq — > 00 that is, from Eq. (|108|) , Po nt (R, R) — > 0. Indeed if in the absence 



of the cavity the two of photon are emitted one toward right and one toward left, 
it is impossible to detect two from the same side. When R > 0.8, 1Z is practically 
independent from R while it presents an oscillation of period tt in x. Observing Fig. 
6a is evident that when R > 0.5, near the resonance x ~ ir we have TZq ~ 2, while 
near x = (2n + l)vr/2, n integer, there is a region for which TZ < 1. This loss of 
coincidence is due to the fact that only the first and the last pairs of diagrams in Fig. 
5. contribute to TZq but in the first, which gives the probability amplitude of observ- 
ing the coincidence, the two amplitudes Crl and CjnRe 2lu}l ^ c interfere destructively 
for x ~ (2n + 1)tt/2 and R ~ 1 causing zero total amplitude. The second of Eqs. 
( |114| ) is shown in Fig. 7a for < 7Z < 1; the plane part of the graph corresponds 
to the bigger than 1 values. First of all we observe an obvious fact: for R = 0, we 
have 7^-00 = 0, that is if the two photons are emitted both on the same way, it is 
impossible to observe a coincidence in absence of a cavity that mixes the directions. 
From the graph it is evident as when R increment, the reflections increment too, 
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and the coincidence probability arise from the zero value. In Fig. 7b the behaviour 
of is shown for IZoo < 50. We note that when x = (2n + l)7r/2 and R — > 1, 
we have IZoo — > oo, that is the probability of observing a pair of photons from a 
side of the cavity, goes to zero. Indeed the second and the third diagram of Fig. 5 
contribute to T^oo but the third diagram, giving the probability of observation two 
parallel photons, goes to zero for x = (2n+ l)7r/2 and R ~ 1 because of destructive 
interference between Crr and Ci J i l Re lwl l c . Therefore for z>l and realistic reflec- 
tivity, the probability of observing coincidence is always bigger then probability of 
observing two photon on the same side of the cavity. 

On the other hand, we have also seen that for z — ► 0, 7^out can be less then 1 for 
R — > 1 and this is certainly the most interesting case to investigate. In Fig. 8 we 
have shown the behaviour of Eq. (|112|) , as function of x and y, for several values z 
and R = .999. For R > 0.8 the dependence from R is negligible. In Fig. 9 we shown 
the contour plot of IZout, between and 1. From this figure it is evident that the 
dark zones, corresponding to < 1 ; have an extension gradually decreasing for z 
increasing, until they disappear for z > 1 (not shown in Fig. 9). It is interesting to 
note that the probability of observing two photon on one side of the cavity is bigger 
than coincidence probability, when within the cavity the two photons are emitted 
along opposite way. Indeed only third and fourth diagrams in Fig. 5 contribute 
to Po nt (R,R), but while in the fourth diagram the two amplitudes are always in 
phase, in the third diagram the two amplitude Crr and CllR(? w1 ^ c can have op- 
posite phase and interfere destructively for R ~ 1. Then if \Crr\ > \Crl\, that is 
if z > 1, either the third or the fourth diagram, give negligible contributions and 
Pout(R, R) ~ 0. Instead if \Crj_\ > \Crr\ (z < 1) the fourth diagram gives a consis- 
tent contribution to Po u t{R, R) and at the same time the first diagram (proportional 
to Crl but in condition of destructive interference) gives a negligible contribution 
to Pout (R,L). 

Another interesting case is that of the resonance in a broad sense, that is cul/c — ttN, 
with N integer: for N odd there is resonance in a strict sense while for N even there 
is anti-resonance. In this case Eq. (|112|) can be simplified and written as 



4ir2 * 2 - (-l) JV 4Fzcosy + 1 

lZo nt [R, L \R, R) = — , (115) 

z A — (— l) JV 2i' z cosy + r 1 



where F = y2R/ (1 + R). In the Fig. 10a a plot of Eq. ( |115|) is shown as function 



of z and y, for R = .5 and iV odd. The analogue plot for N even can be obtained 
translating the plot by an amount tc along the y axes. It is evident that the ratio 
72-Out is always greater then 1 except for a small region centered around y = tt and 
z = 1/2F which disappear for R — > 1. We can always write Eq. (|115 ) as a ratio 
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between two second degree polynomial in z 



Jz - zl)(z - z u ) 

Ko*t(R, L R, R) = AF 2 ) ±£ =(, 116 

[z — z\)\z — z_) 



where we have defined the roots of the two polynomials as 

1 



_e±* 4 = - Fe ^v. (H7) 



Only for y = tt (mod27r) we can have real positive root: 



2Fz - 1 



y=ir 



K 0nt (R,L\R,R) = — . (118) 



2 



In Fig. 10b a plot of Eq. ( |118| ) is shown for values of z near to 1/2F; from this plot 
we can see in detail the "jump" from the pole to the zero. From Fig. lOd we can 
observe that for R — > 1 the pole in F and the zero in 1/2F tend to the common 
value l/v2 compensating each other, so that T^out — 2. The distance between the 
pole and the zero decreases as ~ (1 — R) 2 and already for R = .9 is less than a part 
on a hundred. It is reasonable to think that for higher and more realistic reflectivity, 
it is not possible to generate really a state so well defined to discriminate between 
the pole and the zero. Furthermore the really physical situation is always described 
by a continuous superposition of modes, therefore we think that effective value of 
T^Out(-R, L\R,R) is ~ 2 in all plane y — z. At last we note that when Crr = 
or Crl = 0, we have respectively 



]imKout{R,L\R,R) = ^-^2, (119) 



lim K 0nt {R, L\R, R) = AF 2 ^> 2. (120) 

z— +oo 

From Eqs. ( |119| - |120| ) we deduce that in these conditions is physically indifferent if 
the two photons are emitted in the same or in the opposite way within the cavity. 
We have already obtained the result T^out = 2, when \ip) = \n±) independently 
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from R, corresponding to equal probability of emission of a pair of photon along 
the same way or in opposite way. In the present case we have the same result for 
ul/c = ttN and R ~ 1. This is consistent with the fact that in the limit of total 
reflectivity in which all frequencies satisfy the resonance (in a broad sense) condi- 
tions uj n = nirc/l, n integer, it is impossible to speak of direction (left or right) of 
emission of a photon, because of the two counter-propagating wave that constitute 
a stationary wave within the cavity have exactly the same weight. Finally we note 
that since Po ut (R, R) = Po ut (L,L), from Eqs (PB}P^) follows 



Pout(R,L) — - ■ ^ — . 



2 + TZ 



Out 



Pout(R, R) 



2 + 7?-out 



1211 



Then when TZo nt = 2 we have 



Pout{R,L) 



Pout(R, R) 



1 

4' 



122) 



in qualitative agreement with Ref. |]T0 . 



5.2 Single photon states 

In this section, we started with investigation of two-photon states, because we were 
interested at the leak of symmetry in the photon-number probability distributions. 
Nevertheless the study of one-photon states is not void of interest. Indeed in previous 
subsection we have shown that in discrete mode representation, the interpretation 
of the results, were "contaminated" from factors \/2 generated by mixing between 
normalization of states with one photon for mode (e.g. |1, 1)) and two photons for 
mode (e.g. 1 2, 0) ) . Now we will see that working with one-photon states this mix- 
ing never appear. In the most general form, the one-photon state \4>) generated by 
Inside operators, can be written as 



I duK^F^-ln), 
i=i J 



(123) 



while the analogous state generated by Outside operators, is given by 
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|F a fa);Out) = C2fo)|0) 

dur] a (uj)\F a (uj); Out). 



(124) 



The probability amplitude to find the electromagnetic field, represented by the state 
\4>) within the cavity, in the state \F a (j]); Out), is 



2 r 

(F a fa);Out|0) = / **vl(f*) [^WHL (125) 
i=i J 



The ratio between the probability Po nt (R) of observing a photon behind mirror 2 
and the probability Pout(^) of observing a photon behind mirror 1 is equal to 



n Q nt(R\L) 



Pout(L) 



(Fi(7/);Out| 



(F 2 (77);Out| 



J du rjl 




[^(w) 


+ C L (cu 






J du 772 


MAM 




+ C R {u 


)r 2 (o/ 





(126) 



where we have redefined = Cr(u) and ^(u;) = Cl(u). Exactly as in section 

5.2, if we assume Cr{uj) (Cl(lu)) proportional to the probability amplitude that an 
active medium within the cavity emit a photon of angular frequency u toward right 
(left), each terms into Eq. ( |126j ) admit a clear physical interpretation illustrated in 
Fig. 11. 

Since in discrete mode representation |/i;In) = b\\0), it is evident that passing to 
discrete case each term between square bracket in Eq. ( |126| ), remain formally un- 
changed, without any a/2 factor. 
For sake of consistency now we must impose 



jdu {\C R (cu)\ 2 + \C L (cu)\ 2 } = l. 



(127) 
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C R (u) + C L (u)n(u)e 



+ 



+ 



FIG. 11. Diagrams illustrating the probability amplitudes (reported in the left 
column), relative to Eq. ( |126| ) Here r\(uS) [^O-^)] is the reflection coefficient of mirror 
1 (at the left) [2 (at the right)]. 

In this case the normalized probability Pout(-R) and Pout(£) are given by 

Pout(i?) = TT^mr Pout(L) = i+kLwy (128) 

6 Conclusion 

We have derived some simple relations for electromagnetic field inside and outside 
an optical cavity, using a non-unitary transformation between Inside and Outside 
operators. The convenience of this approach lies in the fact that we do not need to 
know any details of internal processes that generate the two photon, for calculate 
the photon-number probability distribution outside the cavity. Conversely we can 
obtain information on internal processes, by comparing the calculate and measured 
probability distribution. The method is a natural extension to non-unitary trans- 
formation, of the usual method employed, e.g., in the quantum theory of the lossless 
Beam- Splitter. f|§. 
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Appendix 

In the section 5.2 we introduced the three orthonormal states \n±) e \uq) without 
derivation. This will be done in this Appendix. 

We rewritten the commutation relation for operators b\(u) and b 2 (co) in discrete 
mode representation and symmetrical cavity, as 



b 2 {u),b\{uj) 
b 2 (u),b\(") 



Si(w),&tM 



AM, 



(129) 



where we have defined 



A(c) 



1 -R 2 



l-2Rcos{2ul/c)+R 2 ' 



R 



p(u>) = - 2 yi^r cos M/c)- 



(130) 



Using a slightly different notation with respect to section 5.2, we define 



\ n , C; In) = 1|£^|0), n + C = 2, 



'TV. VC\ 



(131) 



where again the factor {1Z\C\)~ 1 I 2 is due to the possible presence of two photon 
on a single mode. Of course, since b\ and b\ does not commute, this kets not form 
an orthonormal base, but they are however linearly independent. Indeed if we define 



{K, C- In|7e', C!\ In) = G(K, C; K', £), 



(132) 
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we can calculate, using Eqs. (|129[ - |130D 



G(TZ,C;TZ',C) = A 2 



/ 1 V2p p 2 \ 

V2p 1 + p 2 V2p 
\P 2 v^P 1 / 



(133) 



Therefore the kets defined in Eq. (|131|) are linearly independent being their Gram 



determinant positive [25]: 



Det[G] = A 6 (l -p 2 ) 3 > 



(134) 



By diagonalization of G, after some algebra we obtain the orthonormal base we look 
for: 



\n+) = ~|2, 0; In) + 1; In) + i|0, 2; In), 



|n ) = -i=|2,0;In) + i=|0,2;In), 



(135) 



|n_> = i|2, 0; In) - ±=\1, 1; In) + ^|0, 2; In). 



Note on matrix G 

The form of the matrix G is particular and justify this little note. Let So, Si and 
S2 the following 3x3 matrix: 



'1 (T 
S = ( 1 

.0 1, 



Si 



'0 f 

1 

1 0, 



'0 1 o x 

1 1 
,0 1 0, 



(136) 
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It can be readily shown that they satisfy the following multiplication table: 



S^ ■ Sj, 


So 


Si 


s 2 


So 


So 


Si 


s 2 


Si 


Si 


So 


s 2 


s 2 


s 2 


s 2 


So + Si 



(137) 



Now consider the generic matrix Af(a) , (a 6 1) given by 

M(a) = S + a 2 Si + V2aS 2 . 

It is characterized by 

Det[JV(a)] = (1 - a 2 ) 3 , Tr[Af(a)\ = 3 + a 2 . 



(138) 



(139) 



If we indicate with A , \± his eigenvalues and with n , n ± the corresponding eigen- 
vectors, we can write 



'n+ = 1(1,^2,1) A+ = (l + a) 
< n = ^(-1,0,1) A = (l-a 2 ; 

k n_ = |(l,-V2,l) A_ = (l-a) 
Using Eq. (|137|) it is easy to see that 

a + (3 



(140) 



Af(a) ■ M{j3) = (l + «/5)W 



1 + a(3 



(141) 



Since M (0) = So is the identity matrix, it is clear that the inverse of N{a) is still 
a matrix of the form ( |138| ). Indeed putting j3 = —a into Eq. (p.41|) we obtain 



M{a) • Af(-a) = (1-« 2 ) 2 S 



(142) 
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that is 



a 



1 — a 



2\2 



Af(-a). 



(143) 



Finally, from Eq. ( |133| ) we get 



G = A 2 AT(p). 



(144) 



As a curiosity, we note that in two dimension the matrices J\f(a) which satisfy the 
algebra ( |141| ) are given by 



J\f(a) = a + aai, 



(145) 



where a = I and <j\ is the first of the Pauli matrices |2^, p. 160]. It is easy to show 
that 



Det[jV(a)] = 1 - a 2 , Tr[M(a)] = 2. 



(146) 



The eigenvalues X± and the corresponding eigenvectors n± are given by 



114 



11 



|(1,1) \ + = l + a, 



(147) 



■1,1) A_ = l 



a. 
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FIGURE CAPTIONS 



FIG. 6. Plots of TZo for different range of values. In (a) the plane part correspond 
to values of 7Z greater than 2.5. 

FIG. 7. Plots of 7^00 for different range of values. In (a) the plane part corre- 
spond to values of TZ^ greater than 1. 

FIG. 8. Four plot of 7^out calculated for a symmetrical cavity, R = 0.999 and 
several values of z. The dependence from R for R > 0.8 is negligible and not 
reported in the figure. 

FIG. 9. Contour plot of corresponding plot in Fig. 8, shows for values 7^out 
between and 1. It is evident that for z increasing, the zone on the plane x — y in 
which 7£out < 1 decreasing. 

FIG. 10. (a) Plot of Eq. (|115|) for R = 0.5. This is not a realistic value, but we 
choose it for show in a clear manner the various quantities, (b) The same that in 
(a) for y = tt; the vertical straight lines passing through the pole in F and through 
the zero in 1/2F. (c) The same that in (a) for z = 1/2F. (d) Plot of 1/2F (up) and 
F (down), as function of the reflectivity R. 
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